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ABSTRACT 

We study filamentary structure in the galaxy distribution at z ^ 0.8 using data from 
the Deep Extragalactic Evolutionary Probe 2 (DEEP2) Redshift Survey and its evo- 
lution to z ~ 0.1 using data from the Sloan Digital Sky Survey (SDSS). We trace 
individual filaments for both surveys using the Smoothed Hessian Major Axis Fila- 
ment Finder, an algorithm which employs the Hessian matrix of the galaxy density 
field to trace the filamentary structures in the distribution of galaxies. We extract 33 
subsamples from the SDSS data with a geometry similar to that of DEEP2. We find 
that the filament length distribution has not significantly changed since z ~ 0.8, as 
predicted in a previous study using a ACDM cosmological N-body simulation. How- 
ever, the filament width distribution, which is sensitive to the non-linear growth of 
structure, broadens and shifts to smaller widths for smoothing length scales of 5 — 10 
h^^ Mpc from z ^ 0.8 to z ^ 0.1, in accord with N-body simulations. 
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1 INTRODUCTION 

The observed large-scale distribution of galaxies shows 
dense linear features: filaments of galaxies which surround 
huge voids that appear largely empt y, while rich clus- 



ters are found at their intersection (iDavis et al 



de Lapparent. Geller. fc Huchralll986l : iBond fc Myerd 
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Gott et al.ll2005l ). These structures are widely believed to 



have evolved through gravitational instability from small 
density fluctuations in the early universe. The evolution of 
large scale structure with cosmic time can probe the com- 
plex physics that governs the creation of galaxies in their 
host dark matter potential wells. A numb e r of studies (e.g. 
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Coil et'al] l2004bl: lOuchi et al.l |2004 iLe Fevre et al.l 
Meneux et al.l 120061 ) have focused on the redshift depen- 
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dence of the galaxy two point correlation function as a crit- 
ical test of both cosmological and galaxy evolution models. 
The comoving correlation length of galaxies is observed to 
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be almost constant with redshift, which is interpreted as a 
consequence of the increasing bias of galaxies with redshift. 

The two-point correlation function is a complete statis- 
tical measure of galaxy clustering only in the linear regime. 
Statistics of galaxy filaments, such as their lengths and 
widths, can be used as another useful tool to measure the 
large scale structure and test both cosmology and galaxy for- 
mation model s. Filaments, with a typi cal length of 50 — 70 
h^^ Mpc (Bharadwai. Bhavsar fc Shet h 2004), have been 
seen in every wide- field redshift survey, from the Great 
Wall of the CF A2 llde Lapparent. Geller. fc HuchTal Il986l ; 
iGeller fc Huchral Il989l ) to the very long filaments found 
("Gott et al."2005l) in the Sloan Digital Sky Survey (SDSS) 
(York et al. 200Q). Features qualitatively similar to the ob- 
served filamentary structures are also seen in numerical cos- 
mological simulations. Various techniques have been pro- 
posed to identify and charact erize filaments in observational 
and simulated samples (e.g. Moody, Turner fc Gott '1983|; 
Eriksen et al] [2004: Stoica et al. 2005 ; Lacostc et al. 2005 



Novikov. Colombi fc Dori I2OO6I; Uragon-Calvo et al.l [200 



Sousbie et all l2008alld; ISoubie. Colombi fc PichonI 



Forero- Romero et al] |2009| ). iBond. Strauss, fc £e3 
2OIOI , hereafter Paper I and II, repectively) use the eigenvec- 
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tors of the Hessian matrix of the smoothed galaxy density 
field to identify and quantify filamentary structures. The 
filament length and width distributions of the observed lo- 
cal galaxy distribution from SDSS are consistent with those 
from N-body simulations at 2 ~ adopting a standard cos- 
mology. The time evolution of the filament network was 
studied in Paper II using cosmological N-body simulations; 
they found that the backbone of the filamentary structure 
is in place at z = 3. These simulations show that non-linear 
growth of structure has little impact on the length of fil- 
amentary structures, but a great deal on the width. The 
dark matter filament width distribution evolves from z ~ 3 
to 2 ~ 0, broadening and peaking at smaller widths as the 
universe expands. 

Although many recent papers have s tudied the filamen- 
tary structure s of local galaxy surveys (|Stoica et al. I l2007l : 
ISousbie et al.l 2008a, b: Gay ct al. 2003), no study has been 
done on the evolution of filaments in redshift survey data. 
This is mainly due to the small volumes and the result- 
ing severe cosmic variance of existing high redshift surveys. 
However, thanks to the successful completion of the Deep 
Extragalactic Evolutionary Pro be 2 (DEEP2) Galaxy Red- 
shift Survey (jPavis et al.ll2003h . we can study the galaxy 
distribution at z ~ 1 over a large comoving volume (5 x 10* 
h~^y['pc?') over four widely separated fields. 

In this study, we present measurements of filament 
statistics both for the galaxy distributi on at z ^ 0.8 us - 
ing the DEEP2 Galaxy Redshift Survey (jPavis et al.ll2003h . 
which is an R-band-selected survey with a sampling density 
comparable to local surveys, and for the local galaxy distri- 
bution from the SDSS. We identify filaments in the galaxy 
distribution in DEEP2 using the Hessian matrix method of 
Paper I and II, and draw subsamples from the SDSS redshift 
survey with the DEEP2 geometry and sampling, to make a 
direct comparison between the two. 

This paper is organized as follows: § [2] summarizes the 
methods we use to find filamentary structures and measure 
their properties, referring the reader to Paper I and Paper 
II for more details. In § [3l we provide details of the data 
samples used here. §|3]presents our results and § [S] discusses 
their meaning and implications. We assume a fiat concor- 
dance ACDM cosmology with Q.m = 0.3, JIa = l~Q.m = 0.7 
and Hq = 100 h km s~^ Mpc~^ throughout this paper. 



2 METHODS 

In this paper, we use the Smoothed Hessian Major Axis Fil- 
ament Finder (SHMAFF), an algorithm that uses the eigen- 
vectors of the Hessian matrix of the smoothed galaxy distri- 
bution to identify filamentary structures of galaxy data. The 
detailed methodology of SHMAFF and its applications are 
described in Papers I and II; Paper I quantifies the promi- 
nence and shapes of structures in the galaxy distribution us- 
ing the Hessian matrix, while Paper II describes a method 
to find individual filaments, and compares their properties 
in cosmological N-body simulations to those in the SDSS 
galaxy distribution. We summarize the basics here. Since the 
geometry of the DEEP2 fields allows us to study filamentary 
structures best in two dimensions, we use a two-dimensional 
version of SHMAFF in this study. 



2.1 Smoothed density field and its Hessian 

To trace individual filaments in the galaxy distribution, we 
generate the density field and its second derivatives; fila- 
ments will be defined as regions with one eigenvalue of the 
Hessian matrix much larger than the other two. The density 
field is smoothed with a Gaussian kernel with smoothing 
length /, 



P(x) 



j /(x-x')p(x')d^x', 



(1) 



where p(x) is the unsmoothed density field, and f{x) is the 
smoothing kernel. In this study we use Gaussian smoothing. 



/(x) = 



-e 



(2) 



where I is the smoothing length and the smoothing is per- 
formed over a two-dimensional box. The unsmoothed den- 
sity field is given as a sum of delta functions at the positions 
of the galaxies and the smoothed Hessian (the matrix of sec- 
ond partial derivatives) is given by 



f(x-x')^^^dV-- 



9V(x-x') , 2 / .ON 

dx'dx'.. ^(")dx.(3) 



Since the filamentary structures can appear on a variety 
of scales, it is important to smooth on a series of length scales 
/ when searching for filaments. 

The Hessian matrix describes the local curvature of the 
density field, that is, the major axis is aligned along the di- 
rection of lowest concavity. We compute the eigenvalues, Ai 
of the Hessian matrix defined such that Ai < A2, and eigen- 
vectors, Ai, which give the orientation of the structure at a 
given grid cell. The direction of lowest concavity is expected 
to be along the filament itself. We thus simply need to find 
the major axis of the Hessian ellipsoid in order to find the 
direction of a filament at a given point in space. 

2.2 Filament-Finding Parameters 

We trace filaments over those grid cells that satisfy the cri- 
teria Al < 0, p > p, where p is the mean density of the 
smoothed density field. We choose a starting point at the 
local maximum density, and trace out the filament both par- 
allel and antiparallel to the A2 axis until its local curvature 
exceeds a given threshold. Along a filament, if the angular 
rate of change of the axis of structure exceeds a value C, 
filament tracing is stopped and the point will be marked as 
a filament end. The stopping condition at pixel m is given 
by 



I > sin(CA), 



(4) 



where A is the distance between pixels. 

As each filament is found, the pixels associated with it 
are removed from further consideration as filament starting 
points. In particular, we introduce another input parameter 
K and define a removal width Wi at grid cell i as 



Al i 



(5) 



All pixels within a removal width Wi of the most recently 
chosen filament element are excluded in the next iteration 
of the filament finding procedure, which prevents filaments 
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from being multiply counted. Here, a filament element is 
defined to be a segment of the filament with length equal to 
the distance between pixels. 

This process depends on three input parameters: the 
smoothing length / {h~^ Mpc), the curvature criterion C 
and the width of filament removal K. For tracing 
filaments in the large-scale galaxy distribution, the best 
values of the input parameters suggested by Paper II are 
C = 30°l~^ and 40°Z~^ on smoothing scales of Z = 5 and 10 

Mpc respectively, and K = 1 for all smoothing scales. 
These parameters were determined in three dimension; we 
adopt these values here in two dimension. 



2.3 Filament Measurement 

The length of each filament is defined as the distance along 
the filament between its two ends, which are specified by 
Eq. |4]or where the density no longer exceeds the threshold. 
As discussed in Paper II, the filament finder would identify 
an isolated spherical over-dense region as a "filament" of 
length of order the smoothing length. Thus in this paper, 
we exclude "filaments" whose lengths are shorter than the 
smoothing length I. 

The width of a filament element (W) is defined to be the 
root mean squared perpendicular offset of galaxies within a 
smoothing length I: 



W = 



yiv |-R,|2 



TV 



(6) 



where the sum is over the A'^ galaxies within a smoothing 
length of the filament element. Here, is the perpendicular 
offset of nearby points from the filament axis, and is defined 
as 

Ri = Aj X (Ai X (xj - Xi)), (7) 
where Aj is the unit vector along the axis of structure. 



3 DATA 

The main data analyzed in this paper come from the DEEP2 
and SDSS surveys. We now describe these two surveys to 
understand the samples and their selection funtions. 




Figure 1. Aitoff projection in equatorial coordinates of the 
ang ular coverage of t lic SDSS Data Release 7 galaxy catalog 
llAbazaiian et al.|[2009l) . Marked with a black solid line is the re- 
gion from which samples used in this paper were drawn (8'' ^ 
a ^ 16'^ and 0° ^ (5 ^ 60°). 



of Rab = 24.10. The survey spans a comoving volume of 
approximately 5 x 10^/i~^Mpc^, covering ~ 3 deg^ over four 
widely separated fields to limit the effect of cosmic variance. 
The DEEP2 observations, catalog const ruction, and d ata re - 



duction are described in more detail in iDavis et all ||2003|) : 
ICoil et al.1 (|2004ah and lOavis. Gerke. &: NewmanI (j20oJ r^ 

In fields 2, 3, and 4, the spectroscopic target galax- 
ies are preselected using a color cut in [B — R)-{R — I) 
space to ensure that most galaxies have redshifts greater 
than 0.75. With this color cut, ~ 90% of the targeted 
galaxies are at z > 0.75, and only ~ 3% of the z > 0.75 
galaxies brighter t han the magnitude limit are not selected 
(jPavis et al.boOSi '). A fourth field, the extended Groth Strip 
(EGS), does not have this redshift preselection. We use the 
absolute B-band magnitu des (Mb) and restfra me U — B 
colors that are derived in IWillmer et al.l l|200(j) ; we apply 
no corrections for luminosity evolution. We create volume- 
limited samples as a function of Mb in three fields covering 
~ 2.2 deg^. The DEEP2 sample is not complete, and has a 
complicated angular mask. Details of sample definitions are 
discussed in Section 13.31 



3.1 DEEP2 Galaxy Redshift Survey 

The high-redshift galaxy sample used i n this paper is from 
the DEEP2 Galaxy Redshift Survey (|Davis et al.l l2003f ). 
which used the Deep Imaging and Multi-Object Spectro- 
graph (DEIMOS) (FabCTct al. 2003) on the 10 m Keck II 
telescope to obtain spectra of optically selected galaxies at 
z ~ 1. The selection was done from deep BRI photome- 
try drawn from images taken with the CFHT12k camera on 
the C anada-France-Hawaii Telescope (CFHT) l|Coil et al.l 
l2004bl ). The DEEP2 spectra have a resolution of i? ~ 5000, 
and rms redshift errors, which are determined from repeated 
observations, are ~ 30 km s~^. The survey has measured 
high-confidence redshifts for ~ 28, 100 galaxies in the red- 
shift range of 0.7 < z < 1.5 down to a limiting magnitude 



3.2 SDSS 

We will compare the DEEP2 filamentary structures to those 
obtained from the SDSS at redshifts oi z < 0.1. The SDSS 
is an extensive photometric and spectroscopic survey, which 
has obtained photometry of a quarter of the sky and spec- 
tra of over 1.6 million o bjects. Imaging is ob t ained in the 
u, p, r, i, and z bands (jFukueita et al.l Il996l : ISmith et al.l 
l2002l : llvezic et al.ll2004l) with a drift-scan camera with 30 
2048x2048 CCDs ("Gunn ct aJ [l998l ) on a dedicated 2.5 
m telescope ijCunn e t al. 200^). Spectra are measured with 
two fiber- fed digital spectrographs on the same telescope. 
Galaxies are selected for spectroscopy based on a magnitude 



^ All magnitudes in DEEP2 data are i n the AB system. For pho- 
tometric details, see lCoil et al. I ll2004bl') . 
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limit jStrauss et al.|[20o3 ). An overview of the data pipelines 
and da ta products is provide d in the Early Data Release 
paper (jStoughton et al.ll2002l l. The galaxy sample used in 
this paper are constru cted from the SDSS Data Release 7 
l|Abazaiian et al.ll2009l . hereafter DR7). As of DR7, the spec- 
tra of ~ 930, 000 galaxies have been measured, covering 9380 
deg'^. Galaxy redshift errors are typically ~ 30km s~^, which 
is similar to DEEP2. 

For this analysis, we make use of the New York Univer- 
sity Value Added Galaxy Catalog (NYU-VAGC), which is a 
compilation of the galaxy catalog from the SDSS DR7, pub- 
licly available at http:/ /sdss. phy sics. nyu.edu/vagc/ A de- 
tailed description can be found in iBlanton et al.li|2005l ). We 
construct volume-limited samples from the northern portion 
(8*^ ^5 Q ^ Kf" and 0° ^ 5 ^ 60°) of the survey as shown 
in Figure [T] We use Mo.i^, the r-band absolute magnitude 
corrected to it s z = 0. 1 value using the K-correction code of 
iBlanton et al.l (l2003bll a nd the luminosity evolution model 
of Blanton et al.l ( 2003d ) to define volume-limited samples. 
We will define 33 subsamples from this SDSS sample volume 
shown in Figure[T]to match the geometry of the DEEP2 sam- 
ples. Details of SDSS sample definitions and redshift range 
are discussed in Section [3.31 





3.3 Galaxy Sample Definitions 

3.3.1 Absolute magnitude cut 

As explained in Paper II, sparse sampling of the galaxy den- 
sity field can impact the filament detection rate. A small 
number of galaxies per smoothing volume can create false 
filament detections and cause real filaments to be missed. 
In order to make the most direct comparison of the DEEP2 
and SDSS galaxy distribution, we need to make sure they 
have the same densities. 

The mean galaxy comoving number density of galax- 
ies with absolute magnitudes below an absolute magnitude 
Mcut is: 

^{M)dM, (8) 

- oo 

where 3>(M) is the luminosit y function o f the galaxy sur- 
vey in question. W e adopt ISchechteil l| 19761 ') fits from 
IBlanton et al] (l2003al) for the SDSS galaxy luminosity func- 
tion and from lFaber et al ] 12007. ') for the DEEP2 luminosity 
function. We adjust the magnitude cut of each survey to 
define galaxy samples with matched number densities. We 
adopt the K-corrected rest-frame absolute magnitude cut 
Mb = -19.47 for DEEP2 and Mo.i^ =-19.7 for SDSS with 
h = 1.0, and for which the two samples have comparable 
mean number densities of nj. ~ 0.008 /i^Mpc"^. These mag- 
nitude cuts are ~ 1 magnitude fainter than the characteristic 
magnitude Al* of each survey. In Figure [2l we show the B- 
band absolute magnitude Mb and redshift of each galaxy in 
the DEEP2 catalog and the regions from which our samples 
were drawn. The DEEP2 sample is volume-limited for blue 
galaxies, but not for red galaxies, due to the selection in 
the obse rved R band, which corresponds to the rest-frame 
UV (see lWillmer et all (|2006l ') for more details on selection 
effects in the sample). Due to this selection effect and the 
lower sampling density beyond z ~ 1, we limit our filamen- 
tary study in this paper to 2 < 0.92. 



Figure 2. (Top) Absolute B-band magnitude (in AB magni- 
tudes, with h = I) vs. redshift of the DEEP2 galaxy catalog. The 
magnitude and redshift limits of the sample used here are shown 
in solid lines (0.75 ^ z ^ 0.92 and Mb ^ -19.47). (Bottom) Ab- 
solute Mo.i^ magnitude (with /i = 1) of SDSS DR7 galaxy catalog 
as a function of redshift. The magnitude and redshift limits of the 
main sample used here are shown as solid lines (0.02 ^ z 0.10 
and ^ —19.7). The limits for the sparser sample used to 

study the effect of the direction of sample extraction ( i|3.3.2|l are 
shown as dotted fines (0.02 ^ 2 ^ 0.15 and Mo.i^ ;£ -20.6). 

Figure [2] shows the absolute magnitude Afo.i^ of each 
galaxy from the SDSS catalog and cuts in magnitude and 
redshift with solid lines . Our SDSS sample consists of 528343 
main sample galaxies (jStrauss et al ] l2002l ) with 0.02 < z < 
0.1 and Mo.ir < -19.7. 

The magnitude cuts were chosen to make the two sam- 
ples have comparable number density n$ based on the lu- 
minosity functions. However, the DEEP2 redshift sample is 
roughly 35% incomplete to its magnitude limit and redshift 
limit 0.75 < z < 0.92, du e to unobserved galaxies and red- 
shift failures jWillmer et al. 2006) (this is in addition to the 
3 % incompleteness of the color selection, as described in 
§[3TT}. We thus apply the DEEP2 angular completeness win- 
dow function to the SDSS sample to have the same number 
density (n = 0.005 h^Mpc'^) as the observed DEEP2 sam- 
ple. We do this by selecting subsamples of SDSS with the 
same geometry as DEEP2, as we describe in Section 13.3.31 
after we describe the DEEP2 geometry. 



3.3.2 Subsample definition 

Each DEEP2 field is much longer in the redshift direction 
than on the sky; the 1 — 2 x 0.5 deg'^ fields used for this 
work span 40 — 80 x 20 Mpc in transverse comoving 
extent, while the range 0.7 < z < 1.0 corresponds to 560 
Mpc comoving in the redshift direction. From this comoving 
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Figure 3. Example of the subsample extraction from tlie DEEP2 
(top) and SDSS (bottom) data. Top : Rcdsiiift-space distribution 
of galaxies in Region 2 of DEEP2, shown as a function of redshift 
and comoving distance along the projected distance across the line 
of sight, assuming a ACDM cosmology. We extract one subsample 
from each region with dimensions of 320 X 40 X 14 (h~^Mpc)^ 
(in red). Bottom : One slice with thickness of 14 Mpc from 
SDSS: extracted subsamples are shown in red. We draw 33 such 
subsamples from the SDSS sample volume. 



volume, we select subsamples with dimension of Li x L2 x L3 
= 320 X 40 X 14 (h^^Mpc)^ =179, 200 /i~^Mpc^ and define 
this as a standard box size. The comoving distance of the 
samples along the line of sight spans 1850 - 2170 Mpc, 
which corresponds to 0.75 < z < 0.92. We extract three 
subsamples, one from each region (see the example in the 
upper panel of Figure O). We exclude the fourth field, the 
extended Groth Strip (EGS), which has a narrower width 
than other fields. 

We can extract a larger number of subsamples with the 
standard box size out of the SDSS sample volume. We divide 
the volume into slices of thickness 14 Mpc, and place 
as many rectangles into each slice as possible. Using this 
method, we obtain 33 subsamples in the redshift range of 
0.02 < z < 0.10. The bottom panel of Figure [3] shows an 
example of subsample extraction in an SDSS slice. 

The long axis of the DEEP2 subsamples is along the 
redshift direction, while it is perpendicular to it in SDSS. 
One might be concerned that redshift distortions due to pe- 
culiar velocities along the line of sight would affect the fila- 
mentary statistics differently in the SDSS and DEEP2 sub- 
samples. We were forced to do this, as the volume-limited 
SDSS sample has an extent in the radial direction less than 
the length of the standard box. We test the effect of the di- 
rection of sample extraction by comparing the filamentary 
properties from the subsamples extracted perpendicular to 
the line of sight to those extracted parallel to the line of 
sight. To do this, we build another sparser sample from the 
SDSS over the redshift range of 0.02 < z < 0.15, with an 
absolute magnitude cut M^o.i < —20.6 (Figure [2|. We ex- 
tracted 135 subsamples perpendicular to the line of sight and 
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Figure 4. The width distributions (top) and the length distribu- 
tion (bottom) of filaments of subsamples perpendicular (in blue) 
and parallel (in red) to the line of sight found with the smoothing 
length I = 5 h-^Mpc, C = 30°^!, and K = 1. 



58 subsamples parallel to the line of sight. The filaments are 
found with the smoothing length I = 5 h~^Mpc, C — 30°/~^, 
and K = 1. We derive the length and width distributions 
for each subsample and calculate the composite distributions 
and its error by calculating the mean and a standard devi- 
ation at each bin. As shown in Figure [4] we found that the 
filament width distributions are essentially indistinguishable 
between the perpendicular and parallel subsamples. This is 
confirmed by the Kolmogorov-Smirnov two-sample test (K-S 
test), which we applied to the two distributions. The length 
distributions of the two are essentially identical, differing by 
only 2 CT by the K-S test. The effect of the direction of the 
subsample extraction on the filament finder and filamentary 
properties are negligible presumably because we smooth the 
density field with a smoothing length that is much larger 
than typical galaxy peculiar velocities. 



3.3.3 Survey completeness 

We have to consider the survey completeness. The DEEP2 
survey spectroscopically targets ~60% of objects that pass 
the apparent magnitude and color cuts mentioned above. 
The redshift success is 73% of those targeted galaxies 
(|Willmer et al.ll2006l ). Therefore, we have successful red- 
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Table 1 . Summary of filamentary property 



Smoothing length 


5 h-^ 


Mpc 


10 h- 


^ Mpc 


data 


DEEP2 


SDSS 


DEEP2 


SDSS 


Total number of filaments 
mean # of filaments per subsample 
mean total length of filaments Mpc) 


97 

32.3 ±4.7 
404.0 ±29.5 


958 
29.0 ±4.9 
376.1 ±53.8 


24 
8.0 ± 1.7 
216.6 ±20.0 


245 
7.42 ± 2.2 
186.5 ±43.4 


width distribution peak /i {h^^ Mpc) 
width distribution cr (h^^ Mpc) 


2.021 ± 0.009 
0.508 ±0.005 


1.924 ±0.006 
0.544 ± 0.004 


3.963 ±0.023 
0.580 ±0.016 


3.802 ± 0.017 
0.808 ± 0.011 



shifts for ~ 50% of all galaxies in the surveyed fields with 
apparent magnitude oi R < 24.1. The sampling rate is a 
complex function of position across each field. In order to 
model this effect, we use the angular window function of the 
DEEP2 survey. The mask gives the completeness at each an- 
gular position; in unobserved regions such as around bright 
stars, the completeness is zero. We project the angular win- 
dow function covering the three DEEP2 regions onto the ge- 
ometry of each box to generate three-dimensional complete- 
ness maps. Then each SDSS subsample with the geometry 
of the DEEP2 standard box is diluted with a completeness 
map randomly selected among the three DEEP2 samples. 

We have further complications because the sampling 
rate is non-uniform due to the necessities of slitmask de- 
sign. Spectra of objects are not allowed to overlap on the 
CCD when observing with multi-object slit masks, there- 
fore, objects that lie near each other in the direction on 
the sky that maps to the wavelength direction on the CCD 
cannot be observed simultaneously. This results in under- 
sampling in the highest density regions on the plane of the 
sky. In order to reduce the effect of this bias, adjoining 
slit masks are positioned approximately a half-mask width 
apart, giving each gala xy at least two chances to be on a 
mask l|Coil et al.l 120081 ). Despite this, the probability that 
a target with nearest neighbor < 10" away is selected for 
spectro scopy is diminished by ^ 25%. Many DEEP2 related 
papers l|Coil et al.ll2004al : IConrov et a l."2005'') model this ef- 
fect by applying the actual DEEP2 mask-making algorithm 
to the mock galaxy catalogs, which throws out some galax- 
ies located close to other galaxies in the sky. This effect on 
the filamentary properties is n egligible, as it is m ore relevant 
on small scale s < 2 h~^Mpc (|Coil et al.ll200d ). In particu- 
lar, ICoil et al.l (|2008| ) have studied the bias in the two-point 
correlation function due to the slitmask effect; they found it 
was 3.5 % at 1 Mpc, 1 % at 5 Mpc (the minimum 
smoothing length we use), and under a percent by 10 
Mpc. 



The final SDSS subsamples have a mean of 992 galaxies 
each, with a standard deviation of 220, comparable to the 
DEEP2 values (1259, 616, and 1095 galaxies in each field). 
The mean number density of galaxies in the DEEP2 and 
SDSS subsamples is n ~ 0.0055 /I'^Mpc"'^ with a standard 
deviation of ~ 0.001 /i^Mpc"^. 



4 RESULTS 

We project the galaxy distribution in each of the three 
DEEP2 subsamples and 33 SDSS subsamples along the 
short axis and measure a two-dimensional density field. We 
smooth each subsample with smoothing lengths I — 5 and 
10 Mpc and run the filament finder on them using 
C = 30°r^ and 40°r^ respectively. The width of filament 
removal is given as if = 1 for all smoothing scales. Figure [S] 
shows the density field smoothed with 1 = 5 Mpc, maps 
of Al and the identified filaments for DEEP2 Region 2. In 
the top two panels, the red bars indicate the direction of the 
axis of structure at each point. The axis of structure aligns 
with the local filamentary structure, and we identify indi- 
vidual filaments using the curvature criterion C = 30° 
and if = 1. In the bottom two panels of Figure [5] we show 
the identified filaments along with the galaxy distribution 
itself in green dots for the two smoothing lengths. 

After excluding filaments shorter than a smoothing 
length we found 97 filaments in the three DEEP2 sub- 
samples, and 958 filaments in the 33 SDSS subsamples, with 
smoothing length of / = 5h~^ Mpc. The mean total length 
of filaments of each subsample is 404.0 Mpc with a stan- 
dard deviation of 29.5 h'^ Mpc for DEEP2, and 376.1 h'^ 
Mpc with a standard deviation of 53.8 Mpc for SDSS. 
This gives a total filament length per unit area of 3.16 x 10~^ 
h Mpc"^ for DEEP2 subsamples, and 2.94 x 10"^ h Mpc"^ 
for SDSS. These results are in excellent agreement; the over- 
all length of filaments in the two cases is indistinguishable. 
We summarize the filamentary properties of DEEP2 and 
SDSS at the two smoothing length scales in Table 1. 

In Figure [S] we show the filaments in Region 3 and 4 of 
DEEP2 and two selected subsamples of SDSS. The grayscale 
map of Al is shown with the identified filaments in red and 
the galaxy distribution in green dots. The filament distri- 
butions of DEEP2 and SDSS subsamples look qualitatively 
similar. We count how many galaxies lie within one smooth- 
ing length (l) of a filament, and thus calculate the fraction 
of galaxies within a smoothing length of a filament for two 
surveys. With smoothing length / = 5 h"^Mpc, 81.2 ± 1.9% 
of DEEP2 galaxies and 82.1 ± 0.5% for SDSS galaxies lie 
within one smoothing length of filaments; the fraction of 
filament galaxies is almost identical in the two surveys. 

Figure [7] shows the length distribution of filament for 
SDSS and DEEP2 with I = 5 Mpc smoothing. There 
are too few filaments at 10 Mpc smoothing to make 
a useful comparison. We derive the length distribution for 
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Figure 5. The filament finder algorithm is shown in action for the subsample that corresponds to DEEP2 Region 2, a 320 X 40 X 14 
(h~^Mpc)^ slice in the redshift range 0.75 < z < 0.92. The smoothing length is I = 5 h~^Mpc. (a) An inverted grayscale density map 
(lighter=more dense) along with a 20% random sample of bars that indicate the direction of the axis of structure at each point, (b) 
a grayscale map of Ai {lighter=more negative), the first eigenvalue of the second partial derivative of the density field, (c) The galaxy 
distribution (green dots) and filaments (in red) are shown with the grayscale map of Ai (lighter=more negative). The filaments are found 
with parameters C = 30°r^ and K = 1.0. (d) Same as (c) but for the smoothing length / = 10 h ^Mpc. The filaments are found with 
parameters C = 40°i~^ and K = 1.0. 



each subsample and calculate the composite distribution and 
its error by calculating the mean and a standard deviation 
at each bin. The length distributions found in DEEP2 are 
si milar to thos e found in SDSS. The numerical simulations 
m iBondl i 20081 ) showed that dark matter filaments have an 
exponential length distribution at large filament lengths that 
very closely matched that found in a Gaussian random field 
with the same power spectrum. This means that even if 
the filaments in the two galaxy distributions are at different 
stages of their evolution, the length distribution should be 
similar between the two. 

However, the width distribution of filament elements 
changes significantly as non-linear evolution proceeds (Pa- 
per II): it broadens and peaks at smaller widths with cos- 
mic time, according to N-body simulations. We show the 
width distributions of filament elements for DEEP2 and 
SDSS in Figure[8]for the two smoothing lengths. We applied 
a K-S two-sample test to compare the two; the calculated 
probability that the two are drawn from the same distri- 
bution is shown in each panel in Figure [S] In both cases. 



the probability is negligibly small. We make Gaussian fits 
to the distributions to compare the two surveys. For I — 5 
Mpc, DEEP2 has a width distribution that peaks at 
fi = 2.021 ± 0.009 Mpc with a = 0.508 ± 0.005/i"^ 
Mpc, and SDSS has ^ = 1.924 ± 0.006 h'^ Mpc with 
a = 0.544 ± 0.004 ft"^ Mpc. In case oi I = 10 h'^ Mpc, 
the distribution has /x = 3.963 ± 0.0023 Mpc with 

a = 0.580±0.016 Mpc for DEEP2, and fi = 3.802±0.017 
h''^ Mpc and a = 0.808 ± 0.011 /i"^ Mpc for SDSS. For 
both smoothing lengths, the filament element width distri- 
butions broaden and shift to smaller widths from z ~ 0.8 
to z ~ 0.1. Figure [8] shows better match between the width 
distributions of SDSS and DEEP2 at high wi dths than low , 
in agreement with the simulations shown by iBondl (|2008l ). 
Thus, the widest filaments narrow more slowly than do the 
narrowest on es, as one would expect in the ellipsoidal col- 
lapse model (|ZerDovichlll970l ). in which overdensities se- 
quentially contract along their principal axes, in order of 
increasing length. 
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Figure 6. The distribution of filaments in the DEEP2 (upper two panels) and two randomly chosen SDSS subsamples (lower two panels). 
The smoothing length is Z = 5 h~^Mpc and the filaments are found with filament-finding parameters C = 30°l~^ and K = 1.0. The 
grayscale map of Ai (lighter=morc negative) is shown along with the galaxy distribution (green dots) and filaments (in red). 
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Figure 7. The composite length distributions of the filaments of 
three DEEP2 subsamples (in blue) and 33 SDSS subsamples (in 
red). 



5 SUMMARY 

We study the time evolution of the filament network in the 
galaxy distribution by comparing the filamentary structure 
at z ~ 0.8 from the Deep Extragalactic Evolutionary Probe 
2 (DEEP2) Redshift Survey and those at z ~ O.f from 
the Sloan Digital Sky Survey (SDSS). We trace individual 
filaments for both surveys using SHMAFF, an algorithm 
which employs the Hessian matrix to trace the filamentary 
structures in the distribution of structure. We define three 
subsamples from DEEP2 and 33 subsamples from SDSS, 
with the same sampling and geometry, namely a box of 
320 X 40 X 14 (h"^Mpc)^. We smooth the galaxy distribu- 
tion with length scales of / = 5 and 10 Mpc, and trace 
individual filaments along the axis of structure, and mark 
the end of filaments when the axis orientation changes more 
rapidly than a preset threshold of C = 30 and 40°/~^ respec- 
tively. We found 97 filaments in DEEP2 subsamples and 957 
filaments in SDSS subsamples with smoothing length 1 — 5 
Mpc, and 24 filaments for DEEP2 and 230 for SDSS for 
I — 10/i~^ Mpc. Thus the number of filaments per unit vol- 
ume is unchanged from high to low redshift. We find that 
filament length distribution has not changed significantly 
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Figure 8. The composite widtli distributions of the filaments of 
three DEEP2 subsamples (in blue) and 33 SDSS subsamples (in 
red). 



since 2 ~ 1, however, the filament width distribution, which 
is sensitive to non-linear growth of structure, broadens and 
shifts to smaller widths for smoothing length scales of 5 and 
10 ft"^ Mpc from z ~ 0.8 to z ~ 0.1. The evolution in the 
length and width distributions is consistent with predictions 
from a ACDM cosmological N-body simulation. As found in 
Paper II, non-linear growth of structure has a great impact 
on the width of filamentary structures. 

We restricted our study to two-dimensional analysis due 
to the the geometry of the DEEP2 survey. In order to bet- 
ter show the filamentary evolution, however, the filamen- 
tary structures should be studied in larger volumes and an- 
alyzed in three dimensions. The next generation of galaxy 
surveys, which will target the early universe with larger 
volume and depth can open up the possibility of detailed 
study of filamentary structures and their evolution. These 
surveys include the Advanced Dark Energy Physics Tele- 
scope (ADEPT), a space-based spectroscopic survey that 
promises to determine the locatio n of 100 million ga lax- 
ies at 1 < 2 < 2, the BigBOSS (|Schlegel et al.ll2009t ). a 
proposed ground-based wide field spectroscopic survey at 
0.2 < z < 3.5, and all hemisphere H I redshift surveys wit h 
the Square Kilometer Array (SKA) (|Rawlings et al.ll2003 ). 
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